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We present a discrete form of the Wheeler-DeWitt equation for quantum gravitation, based 
on the lattice formulation due to Regge. In this setup the infinite-dimensional manifold of 3- 
geometries is replaced by a space of three-dimensional piecewise linear spaces, with the solutions 
to the lattice equations providing a suitable approximation to the continuum wave functional. The 
equations incorporate a set of constraints on the quantum wavefunctional, arising from the triangle 
inequalities and their higher dimensional analogs. The character of the solutions is discussed in the 
strong coupling (large G) limit, where it is shown that the wavefunctional only depends on geometric 
quantities, such as areas and volumes. An explicit form, determined from the discrete wave equation 
supplemented by suitable regularity conditions, shows peaks corresponding to integer multiples of 
a fundamental unit of volume. An application of the variational method using correlated product 
wavefunctions suggests a relationship between quantum gravity in n -|- 1 dimensions, and averages 
computed in the Euclidean path integral formulation in n dimensions. The proposed discrete 
equations could provide a useful, and complementary, computational alternative to the Euclidean 
lattice path integral approach to quantum gravity. 

^e-mail address : HHamber@uci.edu 

^e-mail address : R.M.Williams@damtp.cam.ac.uk 



1 Introduction 



In this paper we will present a lattice version of the Wheeler-DeWitt equation of quantum gravity. 
The approach used here will be rooted in the canonical formulation of quantum gravity, and can 
therefore be regarded as complementary to the Euclidean lattice version of the same theory dis- 
cussed elsewhere. In the following we will derive a discrete form of the Wheeler-DeWitt equation 
for pure gravity, based on the simplicial lattice formulation of gravity developed by Regge and 
Wheeler. It is expected that the resulting lattice equation will reproduce the original continuum 
equation in some suitable small lattice spacing limit. In this formulation the infinite-dimensional 
manifold of 3-geometries is replaced by the space of three-dimensional piecewise linear spaces, with 
solutions to the lattice equations then providing a suitable approximation to the continuum gravi- 
tational wavefunctional. The lattice equations will provide a new set of constraints on the quantum 
wavefunctional, which arise because of the imposition of the triangle inequalities and their higher 
dimensional analogs. The equations arc explicit enough to allow a number of potentially useful 
practical calculations in the quantum theory of gravity, such as the strong coupling expansion, the 
weak field expansion, mean field theory, and the variational method. In this work we will provide 
a number of sample calculations to illustrate the workings of the lattice theory, and what in our 
opinion is the likely physical interpretation of the results. 

In the strong coupling (large G) limit we will show that the wavefunctional depends on geometric 
quantities only, such as areas, volumes and curvatures, and that in this limit the correlation length 
is finite in units of the lattice spacing. An explicit form of the wavefunctional, determined from 
the discrete equation supplemented by suitable regularity conditions, shows peaks corresponding 
to integer multiples of a fundamental unit of volume. Later the variational method will be in- 
troduced, based here on correlated product (Jastrow-Slater type) wavefunctions. This approach 
brings out a relationship between ground state properties of quantum gravity in n + 1 dimensions, 
and certain averages computed in the Euclidean path integral formulation in n dimensions, i.e. 
in one dimension less. Because of its reliance on a diflFerent set of approximation methods, the 
3-1-1 lattice formulation presented here could provide a useful, and complementary, computational 
alternative to the Euclidean lattice path integral approach to quantum gravity in four dimensions. 
The equations are explicit enough that numerical solutions should be achievable in a number of 
simple cases, such as a toroidal regular lattice with N vertices in 3-1-1 dimensions. 

An outline of the paper is as follows. In Section 2, as a background to the rest of the paper, we 
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describe the formalism of classical gravity, as set up by Arnowitt, Dcscr and Misner. In Section 3, 
we introduce the continuum form of the Wheeler-DeWitt equation and, in Section 4, describe how 
it can be solved in the minisuperspace approximation. Section 5 is the central core of the paper, 
where we transcribe the Wheeler-DeWitt equation to the lattice. Practical details of the lattice 
version are given in Section 6 and the equation solved in the strong coupling limit in both 2+1 
and 3+1 dimensions. A general solution at the full range of couplings requires the inclusion of the 
curvature term, which was neglected in the strong coupling expansion, and Sections 7 and 8 outline 
methods of including this term, by perturbation theory and by the variational method. Section 
9 gives a short outline of the lattice weak field expansion as it applies to the Wheeler-DeWitt 
equation. Section 10 concludes with a discussion. 



2 Arnowitt-Deser-Misner (ADM) Formalism and Hamiltonian 

Since this paper involves the canonical quantization of gravity [1], we begin with a discussion of the 
classical canonical formalism derived by Arnowitt, Dcscr and Misner [2]. While many of the results 
presented in this section are rather well known, it will be useful, in view of later applications, to 
recall the main results and formulas and provide suitable references for expressions used later in 
the paper. 

The first step in developing a canonical formulation for gravity is to introduce a time-slicing 
of space-time, by introducing a sequence of spacelike hyper surf aces labeled by a continuous time 
coordinate t. The invariant distance is then written as 

^ _(^^2 ^ g^^dxi'dx" = gijdx'dx^ + 2gijNUxHt - {N^ - gijN^N^)df , (1) 

where {i = 1,2,3) are coordinates on a three-dimensional manifold and r is the proper time, in 
units with c = 1. The relationship between the quantities dr, dt, dxi, N and Ni basically expresses 
the Lorentzian version of Pythagoras' theorem. 

Indices are raised and lowered with gij{x.) = 1, 2, 3), which denotes the three-metric on the 
given spacelike hypersurface, and Ar(x) and Ar'(x) are the lapse and shift functions, respectively. 
These last two quantities describe the lapse of proper time (AT) between two infinitesimally close 
hypersurfaces, and the corresponding shift in spatial coordinate (AT*). It is customary to mark four- 
dimensional quantities by the prefix ^, so that all un- marked quantities will refer to three dimensions 
(and are occasionally marked explicitly by a ^). In terms of the original four-dimensional metric 
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(2) 
(3) 



gf^y one has 
and for its inverse 

/4g00 4gOj^ ^ /_l/iv2 N^/N^ 

which then gives for the spatial metric and the lapse and shift functions 

g,, = % N = (-^5°°)"'^' Ni = Soi . (4) 
For the volume element one has 

4^9 = N^, (5) 

where the latter involves the determinant of the thrcc-mctric, g = det gij. As usual g^^ denotes the 
inverse of the matrix gij. In terms of these quantities, the Einstein-Hilbert Lagrangian of general 
relativity can then be written, up to an overall multiplicative constant, in the following (first-order) 
form 

£ = ^ ^i? = -gijdtiT^^ - NR"" - NiR\ (6) 
up to boundary terms. Here one has defined the following quantities: 

R' = -2VjTr'^ . (7) 

The symbol V* denotes covariant differentiation with respect to the index i using the spatial three- 
metric gij, and ^i? is the scalar curvature computed from this metric. Also note that the affine 
connection coefficients F^j have been eliminated in favor of the spatial derivatives of the metric 
dkgij, and one has defined tt = tt* j. Since the quantities N and do not appear in the tt*^ dtgij 
part, they are interpreted as Lagrange multipliers, and the " Hamiltonian" density 

n = NR^ + NiR^ (8) 

vanishes as a result of the constraints. Varying the first order Lagrangian of Eq. (6) with respect to 
gij, Ni, N and nij one obtains a set of equations which are equivalent to Einstein's field equations. 
First varying with respect to iTij one obtains an equation which can be viewed as defining TTij, 

dtgij = 2iV5-V2 (^j,.. _ y.. + V,- Ni + ViNj . (9) 
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Varying with respect to the spatial metric gij gives the time evolution for TTjj , 

-27V5'i/2(^ifc^i _ i^^ij-^ ^ ^ {V'V^N - g'^ Vfc N) 

+Vfe {n'^N'') - Vfe A^^ n^^ - Vfe tt*^' . (10) 

Finally varying with respect to the lapse (N) and shift (A^*) functions gives 

R'^{9ij,T^ij) = R\9ij,TTij) = , (11) 

which can be viewed as the four constraint equations = ^R^ — ^R = 0, expressed for this 
choice of metric decomposition [1]. The above constraints can therefore be considered as analogous 
to Gauss's law di = V • E = in electromagnetism. 

Some of the quantities introduced above (such as ^R) describe intrinsic properties of the space- 
like hypersurface, while some others can be related to the extrinsic properties of such a hypersurface 
when embedded in four-dimensional space. If spacetime is sliced up (foliated) by a one-parameter 
family of spacclikc hypersurfaces = x^(x*,t), then one has for the intrinsic metric within the 
spacelike hypersurface 

9ij = 9,..X^X^ with Xf = dix^, (12) 
while the extrinsic curvature is given in terms of the unit normals to the spacelike surface, U^, 

Kijix^t) = -{V^U,)X-X'; . (13) 

In this language, the lapse and shift functions appear in the expression 

dtx^" = NU^" + N'Xl" . (14) 

In the following K = g^^Kij = K'^i the trace of the matrix K. 

Now, in the canonical formalism, the momentum can be expressed in terms of the extrinsic 
curvature 

tt'^ = -^{K'^-Kg'^) . (15) 
It is then convenient to define the quantities H and Hi as (here k = SttG) 

H = 2Kg-^'^Uj'IT'^ - ^TT^) - -L^3^ 

Hi = -2Vj4- (16) 
The last two statements are essentially equivalent to the definitions in Eq. (7). 
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In this notation, the Einstein field equations in the absence of sources are equivalent to the 
initial value constraint 

H{x) = Hi{x) = , (17) 
supplemented by the canonical evolution equations for tt*-' and gij. The quantity 

H = j d^x [N{x)H{x) + N\x)Hi{x)^^ (18) 

should then be regarded as the Hamiltonian for classical general relativity. 
When matter is added to the Einstein-Hilbert Lagrangian, 

I[gA] = / d^x^^R{g^,{x)) + , (19) 

where 4>{x) are some matter fields, the action within the ADM parametrization of the metric 
coordinates needs to be modified to 

I[g, TT, ^,iT^,N] = I dt d^x tt'^ dtgij + 'K^dt^-NT- N'T^ . (20) 

One still has the same definitions as before for the (Lagrange multiplier) lapse and shift function, 
namely N = (— %°°)~^/^ and = g^^ '^goj- The gravitational constraints are modified as well, 
since now one defines 

Ti = m^mgijyn + Hf{gijy^,cl>,7r^), (21) 
with the first part describing the gravitational part given earlier in Eq. (16), 

H{gijy^) = Gij,M tt'^'^' - ^^R + 2\^ 

Hi{gijy^) = -25^,•Vfc7^^■^ (22) 

here conveniently re-written using the (inverse of) the DeWitt supermetric 

Gij^ki = ^9~^^'^ {9ik9ji + 9ii9jk + o:gijgki) , (23) 

with parameter a = —1. Note that in the previous expression a cosmological term (proportional 
to A) has been added as well, for future reference. For the matter part one has 

H'^{9ij,T^'^ ,4>,T^cl>) = y/gToo {gij,TT'\(l),T:^) 

Hf{gijy^,cj>,7r^) = -V5^o^(5^i,7r'^</),7r^). (24) 
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We note here that the (inverse of the) DeWitt supermetric in Eq. (23) is also customarily used 
to define a distance in the space of three-metrics gij{x). Consider an infinitesimal displacement of 
such a metric gij — )■ gij + 6gij , and define the natural metric G on such deformations by considering 
a distance in function space 

\\5gf = j d^xN{x)G'^-''\x)5g,,{x)6gki{x) . (25) 

Here the lapse N{x) is an essentially arbitrary but positive function, to be taken equal to one in 
the following. The quantity G'^^'^\x) is the three-dimensional version of the DeWitt supermetric, 

G'^'"' = 1^/5 (5V+5V' + «5V0 , (26) 

with the parameter a of Eq. (23) related to a in Eq. (26) by a = — 2a/(2 + 3a), so that a = — 1 
gives a = — 2 (note that this is dimension-dependent). 



3 Wheeler-DeWitt Equation 

Within the framework of the previous construction, a transition from a classical to a quantum 
description of gravity is obtained by promoting gij, tt'-', H and Hi to quantum operators, with 
gij and ^'-^ satisfying canonical commutation relations. In particular the classical constraints now 
select a physical vacuum state l^*), such that in the source free case 

H\^) =0 Hi\^) = , (27) 

and in the presence of sources more generally 

f I*) =0 fi\^>) = . (28) 

As in ordinary nonrelativistic quantum mechanics, one can choose different representations for the 
canonically conjugate operators g^j and In the functional position representation one sets 

gij{x) ^ gij{x) 7r'^(x) ^ -ih ■ IQttG ■ ——- . (29) 

In this picture the quantum states become wave functionals of the three- metric gij (x) , 

I*) ^ *[5^,(x)]. (30) 

The two quantum constraint equations in Eq. (28) then become the Wheeler-DeWitt equation 
[3, 4, 5] 



with the inverse supermetric given in 3+1 dimensions by 

Gij^kl = \ 9~^^'^ {9ik9jl + 9il9jk — 9ij9kl) , (32) 
and the difFeomorphism (or momentum) constraint 

l^2igij Vfe ^ + ^'[5^i(x)] = . (33) 

This last constraint imphes that the gradient of ^ on the supcrspace of gij^s and (/)'s is zero 
along those directions that correspond to gauge transformations, i.e. diffeomorphisms on the three 
dimensional manifold, whose points arc labeled by the coordinates x. The lack of covariance of 
the ADM approach has not gone away, and is therefore still part of the present formalism. Also 
note that the DeWitt supermetric is not positive definite, which implies that some derivatives with 
respect to the metric have the "wrong" sign. It is understood that these directions correspond to 
the conformal mode. 

A number of basic issues need to be addressed before one can gain a full and consistent under- 
standing of the dynamical content of the theory [6, 7, 8, 9, 10]. These include possible problems of 
operator ordering, and the specification of a suitable Hilbert space, which entails at some point a 
choice for the inner product of wave functionals, for example in the Schrodinger form 

(*|$) = ld,x[g]^*[g,j]^gij] (34) 

where dfi[g] is some appropriate measure over the three-metric g. Note also that the continuum 
Wheeler-DeWitt equation contains, in the kinetic term, products of functional differential opera- 
tors which are evaluated at the same spatial point x. One would expect that such terms could 
produce S^^") (O)-type singularities when acting on the wave functional, which would then have to be 
regularized in some way. The lattice cutoff discussed below is one way to provide such an explicit 
regularization. 

A peculiar property of the Wheeler-DeWitt equation, which distinguishes it from the usual 
Schrodinger equation = ihdt^, is the absence of an explicit time coordinate. As a result 
the r.h.s. term of the Schrodinger equation is here entirely absent. The reason is of course dif- 
feomorphism invariance of the underlying theory, which expresses now the fundamental quantum 
equations in terms of fields gij, and not coordinates. Consequently the Wheeler-DeWitt equation 
contains no explicit time evolution parameter. Nevertheless in some cases it seems possible to as- 
sign the interpretation of "time coordinate" to some specific variable entering the Wheeler-DeWitt 
equation, such as the overall spatial volume or the magnitude of some scalar field [9] . 
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We shall not discuss here the connection between the Wheeler-DeWitt equation and the Feyn- 
man path integral for gravity. In principle any solution of the Wheeler-DeWitt equation corresponds 
to a possible quantum state of the universe. A similar situation already arises, of course in much 
simpler form, in nonrelativistic quantum mechanics [11]. The effects of the boundary conditions 
on the wavefunction will then act to restrict the class of possible solutions; in ordinary quantum 
mechanics these are determined by the physical context of the problem and some set of external 
conditions. In the case of the universe the situation is far less clear, and in many approaches some 
suitable set of boundary conditions need to be postulated, based on general arguments involving 
simplicity or economy. One proposal [12] is to restrict the quantum state of the universe by requir- 
ing that the wave function ^ be determined by a path integral over compact Euclidean metrics. 
The wave function would then be given by 

^9ijA] = j[dgf,u][d(t>] exp(^-i[g^^,(f>]) , (35) 

where / is the Euclidean action for gravity plus matter 

I = J d^x^{R-2X) - -^J d^x^K - I d^x^Cm . (36) 

The semi-classical functional integral would then be restricted to those four-metrics which have 
the induced metric gij and the matter field (j) as given on the boundary surface S. One would 
then expect (as in the case in nonrelativistic quantum mechanics, where the path integral with 
a boundary surface satisfies the Schrodingcr equation), that the wavefunction constructed in this 
way would also automatically satisfy the Wheeler-DeWitt equation, and this is indeed the case. 



4 Minisuperspace 

Minisuperspace models can in part provide an additional motivation for our later work. The quan- 
tum state of a gravitational system is described, in the Wheeler-DeWitt framework just introduced, 
by a wave function ^ which is a functional of the three-metric gij and the matter fields 0. In general 
the latter could contain fields of arbitrary spins, but here we will consider for simplicity just one 
single component scalar field ^{x). The wavefunction will then obey the zero energy Schrodinger- 
like equation of Eqs. (31) and (33). The quantum state described by * is then a functional on the 
infinite-dimensional manifold W consisting of all positive definite metrics gij{x) and matter fields 
4>{x) on a spacelike three-surface S. We note here that on this space there is a natural metric 

ds''[5g,5ct>\ = I [G^^'^'^ix^x') 5gij5gki{x') + ^5\x-x') 5ct>{x)5(t>{x')\ , (37) 



where 

G'^'''^{x,x') = G'^'''^{x)6^{x-x') (38) 



and 



= i V5 g'^Wix) + g''{x)g^'{x) - 2g'\x)g'\x) (39) 



is the DeWitt supermetric. 

In general the wavefunction for all the dynamical variables of the gravitational field in the 
universe is difficult to calculate, since an infinite number of degrees of freedom are involved: the 
infinitely many values of the metric at all spacetime points, and the infinitely many values of the 
matter field at the same points. One option is to restrict the choice of variable to a finite number 
of suitable degrees of freedom [13, 14, 15, 16, 17]. As a result the overall quantum fluctuations are 
severely restricted, since these are now only allowed to be nonzero along the surviving dynamical 
directions. If the truncation is severe enough, the transverse-traceless nature of the graviton fluc- 
tuation is lost as well. Also, since one is not expanding the quantum solution in a small parameter, 
it can be difficult to estimate corrections. 

In a cosmological context, it seems natural to consider initially a homogeneous and isotropic 
model, and restrict the function space to two variables, the scale factor a{t) and a minimally coupled 
homogeneous scalar field ^(i) [16]. The space-time metric is given by 

dr^ = N^{t)df - gijdx'dx^ . (40) 

The three-metric gij is then determined entirely by the scale factor a{t), 

gij = a^{t) gij (41) 

with gij a time-independent reference three-metric with constant curvature, 

^Rijki = k (gij gki - gu gjk) , (42) 

and A; = 0, ±1 corresponding to the flat, closed and open universe case respectively. In this case 
the minisuperspace W is two-dimensional, with coordinates a and ^, and supermetric 

ds^[a,(j)\ = N-^{-ada^ + dcf)'^) . (43) 

From the above expression for (is^[a, 0] one obtains the Laplacian in the above metric, required for 
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the kinetic term in the Wheeler-DeWitt equation, ^ 

4V>,^) = ^|^a---^| . (44) 

Since the space is homogeneous, the diffeomorphism constraint is triviahy satisfied. Also N is 
independent of gij so in the homogeneous case it can be taken to be a constant, conveniently 
chosen as AT = 1. 

It should be clear that in general the quantum behavior of the solutions is expected to be quite 
different from the classical one. In the latter case one imposes some initial conditions on the scale 
factor at some time t^, which then determines a{t) at all later times. In the minisuperspacc view 
of quantum cosmology one has to instead impose a condition on the wavepacket ^ at a = 0. Due 
to their simplicity, in general it is possible to analyze the solutions to the minisuperspacc Wheeler- 
DeWitt equation in a rather complete way, given some sensible assumptions on how *(a, (j)) should 
behave, for example, when the scale factor a approaches zero. 

In concluding the discussion on minisuperspacc models as a tool for studying the physical 
content of the Wheeler-DeWitt equation it seems legitimate though to ask the following question: 
to what extent can results for these very simple models which involve such a drastic truncation 
of physical degrees of freedom, be ultimately representative of, and physically relevant to, what 
might, or might not, happen in the full quantum theory? 



5 Lattice Hamiltonian for Quantum Gravity 

In constructing a discrete Hamiltonian for gravity one has to decide first what degrees of freedom one 
should retain on the lattice. There are a number of possibilities, depending on which continuum 

theory one chooses to discretize, and at what stage. So, for example, one could start with a 
discretized version of Cartan's formulation, and define vicrbeins and spin connections on a flat 
hypercubic lattice. Later one could define the transfer matrix for such a theory, and construct a 
suitable Hamiltonian. 

Another possibility, which is the one we choose to pursue here, is to use the more economical 

(and geometric) Regge- Wheeler lattice discretization for gravity [18, 19], with edge lengths suitably 

defined on a random lattice as the primary dynamical variables. Even in this specific case several 

^The ambiguity regarding the operator ordering of p^/o = a'^'^'^^^pa'^p in the Wheeler-DeWitt equation can in 
principle be retained by writing for the above operator V'^ the expression 

but this does not seem to affect the qualitative nature of the solutions. The case discussed in the text corresponds 
to g = 1, but g = seems even simpler. 
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avenues for discretization are possible. One could discrctize the theory from the very beginning, 
while it is still formulated in terms of an action, and introduce for it a lapse and a shift function, 
extrinsic and intrinsic discrete curvatures etc. Alternatively one could try to discrctize the contin- 
uum Wheeler-DeWitt equation directly, a procedure that makes sense in the lattice formulation, 
as these equations are still given in terms of geometric objects, for which the Regge theory is very 
well suited. It is the latter approach which we will proceed to outline here. 

The starting point for the following discussion is therefore the Wheeler-DeWitt equation for 
pure gravity in the absence of matter, Eq. (31), 

{ - (""'O)' i,,(.)L,(x) - ( - ) = " 

and the diffeomorphism constraint of Eq. (33) , 

I 2i5,,(x) Vfe(x) I *[5,,(x)] = . (46) 

Note that these equations express a constraint on the state l^*) at every x, each of the form 
H{x.) 1^-) = and Hi (x)]*!-) = 0. 

On a simplicial lattice [20, 21, 22, 23, 24] (see for example [25], and references therein, for a 
more complete discussion of the lattice formulation for gravity) one knows that deformations of the 
squared edge lengths are linearly related to deformations of the induced metric. In a given simplex 
a, take coordinates based at a vertex 0, with axes along the edges from 0. The other vertices are 
each at unit coordinate distance from (see Figures 1,2 and 3 for this labelling of a triangle and 
of a tetrahedron). In terms of these coordinates, the metric within the simplex is given by 

9ij{^) = h (I'oi + loj-llj) ■ (47) 

Note also that in the following discussion only edges and volumes along the spatial direction are 
involved. It follows that one can introduce in a natural way a lattice analog of the DeWitt super- 
metric of Eq. (26), by adhering to the following procedure. First one writes for the supermetric in 
edge length space 

WSl'^f = G'Kl'^) Sl^Sl] , (48) 

ij 

with the quantity G^^{P) suitably defined on the space of squared edge lengths [26, 27]. Through 
a straightforward exercise of varying the squared volume of a given simplex a in d dimensions 

^'(^) = (ir)'det<7,,(i2(a)) (49) 
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to quadratic order in the metric (on the r.h.s.), or in the squared edge lengths belonging to that 
simplex (on the l.h.s.), one finds the identity 

^ E = i.\fd^) [ g'^g'^Sg^^Sgki - g'^g'^jkSgu ] • (50) 

The right hand side of this equation contains precisely the expression appearing in the continuum 
supermetric of Eq. (26) (for a specific choice of the parameter a = —2), while the left hand side 
contains the sought-for lattice supermetric. One is therefore led to the identification 

It should be noted that in spite of the appearance of a sum over simpliccs a, G^^{P) is quite local 
(in correspondence with the continuum, where it is ultra-local), since the derivatives on the right 
hand side vanish when the squared edge lengths in question are not part of the same simplex. The 
sum over a therefore only extends over those few tetrahedra which contain either the i or the j 
edge. 

At this point one is finally ready to write a lattice analog of the Wheeler-DeWitt equation for 
pure gravity, which reads 



-{16^GfGij{l')^f^ - ^) [^R{1^) - 2A] I *[Z2] = , 



(52) 



with Gij{P) the inverse of the matrix G'^^ {t^) given above. The range of the summation over i and 
j and the appropriate expression for the scalar curvature, in this equation, are discussed below and 
made explicit in Eq. (53). 

It should be emphasized that, just like there is one local equation for each spatial point x in the 
continuum, here too there is only one local (or semi-local, since strictly speaking more than one 
lattice vertex is involved) equation that needs to be specified at each simplex, or simplices, with 
Gij defined in accordance with the definition in Eq. (51). On the other hand, and again in close 
analogy with the continuum expression, the wavefunction ^[t^] depends of course collectively on 
all the edge lengths in the lattice. The latter should therefore be regarded as a function of the whole 
simplicial geometry, whatever its nature might be, just like the continuum wavefunction ^[fifij] will 
be a function(al) of all metric variables, or more specifically of the overall geometry of the manifold, 
due to the built-in diffemorphism invariance. On the side we note here that the lattice supermetric 
is dimensionful, Gij ~ Z^"*^ and G^^ ~ Z*^"^ in d spacetime dimensions, so it might be useful and 
convenient from now on to explicitly introduce a lattice spacing a (or a momentum cutoff A = 1/a) 
and express all dimensionful quantities {G,X,li) in terms of this fundamental lattice spacing. 
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As noted, Eqs. (31) or (52) both express a constraint equation at each "point" in space. Here 
we win elaborate a bit more on this point. On the lattice, points in space are replaced by a set 
of edge labels i, with a few edges clustered around each vertex, in a way that depends on the 
dimensionality and the local lattice coordination number. To be more specific, the first term in 
Eq. (52) contains derivatives with respect to edges i and j connected by a matrix element Gij which 
is nonzero only if i and j are close to each other, essentially nearest neighbor. One would therefore 
expect that the first term could be represented by just a sum of edge contributions, all from within 
one {d — l)-simplex a (a tetrahedron in three dimensions). The second term containing ^i2(Z^) in 
Eq. (52) is also local in the edge lengths: it only involves a handful of edge lengths which enter into 
the definition of areas, volumes and angles around the point x, and follows from the fact that the 
local curvature at the original point x is completely determined by the values of the edge lengths 
clustered around i and j. Apart from some geometric factors, it describes, through a deficit angle 
6h, the parallel transport of a vector around an elementary dual lattice loop. It should therefore 
be adequate to represent this second term by a sum over contributions over all {d — 3)-dimensional 
hinges (edges in 3+1 dimensions) h attached to the simplex a, giving therefore in three dimensions 

I -(167rG)2 J2 Gij{a)-^ - 2n,n E + 2A K i = . (53) 

[ i,iCa ^'"J hCa J 

Here Sh is the deficit angle at the hinge h, Ih the corresponding edge length, Va = Vdi'^) volume 
of the simplex (tetrahedron in three spatial dimensions) labeled by a. Gij (a) is obtained either 
from Eq. (51), or from the lattice transcription of Eq. (23) 

Gij,ki{(y) = ^9~^^'^{<y)[9ik{(^)gji{(T) + gii{(j)gjk{cr) - gij{a)gki{a)] , (54) 

with the induced metric gij{cr) within a simplex a given in Eq. (47). The combinatorial factor ng-zi 

ensures the correct normalization for the curvature term, since the latter has to give the lattice 

version of / y^^i? = 2 J2f^ 6hlh (in three spatial dimensions) when summed over all simplices a. The 

inverse of n^^/j counts therefore the number of times the same hinge appears in various neighboring 

simplices, and consequently depends on the specific choice of underlying lattice structure; for a 

fiat lattice of equilateral triangles in two dimensions ricrh = 1/6 . ^ The lattice Wheeler-DeWitt 

equation given in Eq. (53) is the main result of this paper. 

It is in fact quite encouraging that the discrete equation in Eqs. (52) and (53) is very similar to 

what one would derive in Regge lattice gravity by doing the 3+1 split of the lattice metric carefully 

''instead of the combinatorial factor Uah one could insert a ratio of volumes Vcrh/Vfi (where Vn is the volume per 
hinge [23] and Vah is the amount of that volume in the simplex cr), but the above form is simpler. 
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from the very beginning [28, 29, 30]. These authors also derived a lattice Hamiltonian in three 
dimensions, written in terms of lattice momenta conjugate to the edge length variables. In this 
formulation the Hamiltonian constraint equations have the form 

- JLVm^P = 0, (55) 

with Hn defined on the lattice site n. The sum Y^aen extends over neighboring tetrahedra labelled 
by a, whereas the sum J^p&n extends over neighboring edges, here labeled by /3. g[J^ is the inverse 
of the DeWitt supermetric at the site a, and the deficit angle around the edge f3. ^Jg^ is the 
dual (Voronoi) volume associated with the edge /?. 

The lattice Wheeler-DeWitt equation of Eq. (52) has an interesting structure, which is in part 
reminiscent of the Hamiltonian for lattice gauge theories. The first, local kinetic term is the 
gravitational analog of the electric field term E\. It contains momenta which can be considered 
as conjugate to the squared edge length variables. The second local term involving ^i?(/^) is the 
analog of the magnetic (V x ^a)^- In the absence of a cosmological constant term, the first and 
second term have opposite sign, and need to cancel out exactly on physical states in order to give 
i7(x) = 0. On the other hand, the last term proportional to A has no gauge theory analogy, and 
is, as expected, genuinely gravitational. 

It seems important to note here that the squared edge lengths take on only positive values 
> 0, a fact that would seem to imply that the wavefunction has to vanish when the edge lengths 
do, ^{p' = 0) ~ 0. This constraint will tend to select the regular solution close to the origin in 
edge length space, as will be discussed further below. In addition one has some rather complicated 
further constraints on the squared edge lengths, due to the triangle inequalities. These ensure that 
the areas of triangles and the volumes of tetrahedra are always positive. As a result one would 
expect an average soft local upper bound on the squared edge lengths of the type Z| < Zq where Zq 
is an average edge length, (Z?) = l\. The term "soft" refers to the fact that while large values for 
the edge lengths are possible, these should nevertheless enter with a relatively small probability, 
due to the small phase space available in this region. In any case, the nature of the discrete 
Wheeler-DeWitt equation presented here is explicit enough so that these, and other related, issues 
can presumably be answered both satisfactorily and unambiguously. 

The above considerations have some consequences already in the strong coupling limit of the 
theory. For sufficiently strong coupling (large Newton constant G) the first term in Eq. (52) is 



1 



(tr TT^ 



i(tr 7r)2 
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dimension 


dimension of Laplacian 


d dimensions 




2+1 dimensions 


A/l'^ ~ l/A 


3+1 dimensions 


l/l" ~ 1//^ ~ l/V 



Table I: Dimension of the Laplacian term in d dimensions. 



diiiiciisioii 


G dini(nisioii 


A dim(Misi())i 


diuKMisioual 


dim(nisioulcss 


d dimensions 


ld-2 


l/P 


G/VX ~ l'^-^ 


G2/{^-2)A 


2+1 dimensions 


I 


l/P 


G/VX ~ 


G^X 


3+1 dimensions 


P 


l/P 


G/VX - Vt 


GX 



Table 11: Dimensions of couplings in d dimensions. 

dominant, which shows again some similarity with what one finds for non-Abehan gauge theories 
for large coupling g"^. One would then expect both from the constraint h > and the triangle 

inequalities, that the spectrum of this operator is discrete, and that the energy gap, the spacing 
between the lowest eigenvalue and the first excited state, is of the same order as the ultraviolet cut- 
off. Nevertheless one important difference here is that one is not interested in the whole spectrum, 
but instead just in the zero mode. 

Irrespective of its specific form, it is in general possible to simplify the lattice Hamiltonian 
constraint in Eqs. (52) and (53) by using scaling arguments, as one does often in ordinary non- 
relativistic quantum mechanics (for a list of relevant dimensions see Table I and Table II). After 
setting for the scaled cosmological constant A = SttG Aq and dividing the equation out by common 
factors, it can be recast in the slightly simpler form 

[-'"'■ vm°"^'''' ok - = " • 

where one finds it useful to define a dimensionless Newton's constant, as measured in units of the 
cutoff G = IGirG/a?, and a dimensionless cosmological constant Aq = Aqc^, so that in the above 
equation one has a = G/Xq and P = 1/GXq. Furthermore the edge lengths have been rescaled so 
as to be able to set Aq = 1 in lattice units (it is clear from the original gravitational action that 
the cosmological constant Aq simply multiplies the total spacetime volume, thereby just shifting 
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around the overall scale for the problem). Schematically Eq. (56) is therefore of the form 

|-GA, - l^R{s) + l| ns] = , (57) 

with As a discretized form of the covariant super-Laplacian, acting locally on the function space 
of the s = variables. 

We shall not discuss the lattice implementation of the diffeomorphism (or momentum) constraint 
in Eq. (46) . It can be argued that this will be satisfied automatically for a regular or random 
homogeneous lattice. This will indeed be the case for the examples we will be discussing below. 



6 Explicit Setup for the Lattice Wheeler-DeWitt Equation 

In this section, we shall establish our notation and derive the relevant terms in the discrete Wheeler- 
DeWitt equation for a simplex. 

6.1 2+1 dimensions 

The basic simplex in this case is of course a triangle, with vertices and squared edge lengths labelled 
as in Figure 1. We set /qi = ^12 = ^) ^02 = c. 




Figure 1. A triangle with labels. 

The components of the metric for coordinates based at vertex 0, with axes along the 01 and 02 
edges, are 

511 = a, 912 = ^{a + c-b), g22 = c. (58) 
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The area A of the triangle is given by 



= ^[2{ab + bc + ca)-a'^-b'^-c^], (59) 




G'' = A I -1 1 -1 h (60) 



so the supermetric G*-^ , according to Eq. (51), is 

AA ' 

with inverse 

Gij = -2A I 1 1 I . (61) 

Thus for the triangle we have 




^^'^ dsi dsj [dadb db dc dcdaj ^ ^ 



. . / ^2 ^2 52 

Si 5s j 

and the Wheeler-DeWitt equation is 

' + _ 

da db db dc ' dc da 



(16-G)^ ( ^ + ^ + ^1 -2n.,Y.^, + 2Aa| ^s] = 0, (63) 



where the sum is over the three vertices h of the triangle. The combinatorial factor ricrh ensures 
the correct normalization for the curvature term, since the latter has to give the lattice version of 
/ ^ "^R = 2J2h^h when summed over all simplices (triangles in this case) a. The inverse of n^-zj 
counts therefore the number of times the same vertex appears in various neighboring triangles, and 
consequently depends on the specific choice of underlying lattice structure. 
Alternatively, we can evaluate Gij^ki gg.^Qg^^ directly, using 

Gij^ki = -^^{gikgji + gugjk - "^gijgki) (64) 

(note the different coefficient of the last term in two dimensions), with the metric gij as found 
above. The derivatives with respect to the metric are expressed in terms of derivatives with respect 
to squared edge lengths by 

d ^dsm d 



dgijis) "^dgijds. 

This leads to 



d d d 

= ^ + (66) 



^512 5^21 db 
and 



dgn da db 

' ' " (67) 



dg22 db dc 
This procedure gives exactly the same expression for the kinetic term 



d d d 

+ ^- (68) 
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Figure 2. Neighbors of a given triangle. The above picture is supposed to ihustrate the fact that the Laplacian 
A;2 appearing in the kinetic term of the lattice Wheeler-DeWitt equation (here in 2+1 dimensions) contains 
edges a, b, c that belong both to the triangle in question, as well as to several neighboring triangles (here 
three of them) with squared edges denoted sequentially by si = . . . se = /g- 

6.2 3+1 dimensions 

In this case, both methods described for 2+1 dimensions can be followed, but one is much easier 
than the other. 




Figure 3. A tetrahedron with labels. 
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For ease of notation, we define Iqi = a, = b, Iq2 = c, /qs = if^ = e, = /• the 
tetrahedron labehed as in Figure 3, we have 



gn = a, 922 = c, 933 = d, 
912 = ^{a + c - b), gi3 = ^{a + d - e), 523 = ^{c + d - f) 



and its volume V is given by 



(69) 
(70) 



= —[af{-a- f + b + c + d + e) + bd{-b-d + a + c + e + f) + 
+ ce(— c — e + a + b + d + f) — abc — ade — bef — cdf] . 



(71) 



The matrix G*-' is then given by 

/ -2/ e + f-b b + f-e d + f-c c + f-d 

e + f — b — 2e b + e — f d + e — a q 

b + f-e b + e-f -2b r b + c-a 

d + f — c d + e — a r —2d c + d—f 

c + f — d q b + c — a c + d — f —2c 

p a+e—d a+b—c a+d—e a+c—b 



Qt3 



1 

241^ 



p 




a + e — 


d 


a + b- 


c 


a + d — 


e 


a + c — 


b 



\ 



-2a J 



where 



(72) 



p = -2a -2f + b + c + d + e, q = -2c -2e + a + b + d + f, r = -2b -2d + a + c + e + f. (73) 



It is nontrivial to invert this (although it can be done), so instead of using Gij gf.gg . , we evaluate 



with 



G. 



ij,kl 



2V5 



{9ik9ji + 9il9jk — 9ij 9ki) 



d _ A A A 
dgn da db de 

d d d d 

dg22 db dc df 

d _ _^ _^ A 

dg33 dd de df 

d _ _d_ 

dgi2 db 

d _ _d_ 
dgi3 de 

d _ _d_ 
dg23 df 



(74) 



(75) 



The matrix representing the coefficients of the derivatives with respect to the squared edge lengths 
is given in the Appendix, and is the inverse of G*-' found earlier. This is a nontrivial result as it 
acts as confirmation of the Lund-Regge expression which was derived in a completely different way. 
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Then in 3+1 dimensions, the discrete Wheeler-DeWitt equation is 

|-(167rG)2Gi,-^^ - 2n^hY^^Sh + 2Ay|*[s] = 0, (76) 

where the sum is over hinges (edges) h in the tetrahedron. Note the mild nonlocahty of the equation 
in that the curvature term, through the deficit angles, involves edge lengths from neighboring 
tetrahedra. In the continuum, the derivatives also give some mild nonlocality. 

The discrete Wheeler-DeWitt equation is hard to solve analytically, even in 2+1 dimensions, 
because of the complicated dependence on edge lengths in the curvature term, which involves 
arcsin or arccos of convoluted expressions. When the curvature term is negligible, the differential 
operators may be transformed into derivatives with respect to the area (in 2+1 dimensions) or the 
volume (in 3+1 dimensions) and solutions found for the wave function, Figures 4 and 5 give 
a pictorial representation of lattices that can be used for numerical studies of quantum gravity in 
2+1 and 3+1 dimensions, respectively. 




Figure 4. A small section of a suitable spatial lattice for quantum gravity in 2+1 dimensions. 
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Figure 5. A small section of a suitable spatial lattice for quantum gravity in 3+1 dimensions. 
6.3 Solution of the triangle problem in 2+1 dimensions 

In this section we will consider the solution of the Wheeler-DeWitt equation for a single triangle. 
The present calculation is a necessary starting point and should provide a basic stepping stone for 
the strong coupling expansion in 1 jG. In addition it will show the physical nature of the wavefunc- 
tion solution deep in the strong coupling regime. Note that for 1/G ^ the coupling term between 
different simplices, which is due to the curvature term, disappears and one ends up with a com- 
pletely decoupled problem, where the edge lengths in each simplex fluctuate independently. This is 
of course quite analogous to what happens in gauge theories on the lattice at strong coupling, the 
chromo-electric field fluctuates independently on each link, giving rise to short range correlations, 
a mass gap and confinement. Here it is this single-simplex probability amplitude that we will set 
out to compute. 

As in the Euclidean lattice theory of gravity, it will be convenient to factor out an overall 
scale from the problem, and set the (un-scaled) cosmological constant equal to one [23] (see Table 
II). Recall that the Euclidean path integral weight contains a factor Piy) oc exp(— AqF) where 
y = / is the total volume on the lattice. The choice Aq = 1 then fixes this overall scale once 
and for all. Since Aq = 2A/167rG one then has A = SttG in this system of units. In the following 
we will also find it rather convenient to introduce the scaled coupling A 



so that for Aq = 1 (in units of the UV cutoff, or of the fundamental lattice spacing) one has 
A = l/647rG. 




(77) 
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Moreover, it will often turn out to be desirable to avoid large numbers of factors of 167r's by the 
replacement, which we will follow from now on in this section, of 167rG — )■ G. Then A = 1/4G is 
the natural expansion parameter. Note that the coupling A has dimensions of length to the minus 
four, or inverse area squared, in 2+1 dimension, and length to the minus six, or inverse volume 
squared, in 3+1 dimensions. 

Now, from Eq. (63), the Wheeler-DeWitt equation for a single triangle and constant curvature 
density R reads 

r o . . / 52 52 52 



+ 



+ 



+ {2\-R)A^)^[s] 



0, 



(78) 



dh dc dc da 

where a, b, c are the three squared edge lengths for the given triangle, and A/\ is the area of the 
same triangle. In the following we will take for simplicity R = 0. Equivalently one needs to solve 

r 92 92 Q2 



+ 



+ 



da dh dh dc dc da 



+ \\^[a,h,c\ = 



If one sets 

then one can show that 

^2 . 1 



da dh 



* = 



{b + c — a) {a + c — h) 



d2$ 



1 d$ 



+ 



1 



16Aa dAA 



Summing the partial derivatives leads to the equation 

+ 16 A ^ 



dAl 



dAA 



Solutions to the above equation are given by 



^[a,6,c] 



const. — — exp 



±i-4AAy\ 



or alternatively by 



*[a,6,c] = 



^A 



ci cos AAaV X + C2 sin AAaV X 



(79) 



(80) 



(81) 



(82) 



(83) 



(84) 



Note the remarkable, but not entirely unexpected, result that the wavcfunction only depends on 
the area of the triangle Aa(o, 6, c). In other words, it depends on the geometry only. Regularity 
of the wavcfunction as the area of the triangle approaches zero, Aa — >■ 0, requires for the constant 
ci = 0. Therefore the correct quantum-mechanical solution is unambiguously determined, 

1 1 



*[a,5,c] 



V27n7l 



T Aa 



sin 4^aV'^ 



(85) 
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The overall normalization constant has been fixed by the standard rule of quantum mechanics, 

/■oo 

/ dAA I *(^a) P = 1 • (86) 

Jo 

Moreover we note that a bare A < is not possible, and that the oscillatory nature of the wavefunc- 
tion is seen here to give rise to well-defined peaks in the probability distribution for the triangle 
area, located at 

^ (87) 

with n integer. 

6.4 Solution of the tetrahedron problem in 3+1 dimensions 

In this section we will consider the nature of quantum-mechanical solutions for a single tetrahedron. 
Now, from Eq. (76), the Wheeler-DeWitt equation for a single tetrahedron with a constant curvature 
density term R reads 

l^-{16TrGfGij-^^ + {2X-R)V^^s] = 0, (88) 

where now the squared edge lengths si...Sq are all part of the same tetrahedron, and Gij is given 
by a rather complicated, but explicit, 6x6 matrix given earlier. 

As in the 2-1-1 case discussed in the previous section, here too it is found that, when acting 
on functions of the tetrahedron volume, the Laplacian term still returns some other function of 
the volume only, which makes it possible to readily obtain a full solution for the wavefunction. In 
terms of the volume of the tetrahedron Vr one has the equivalent equation for *[s] = /{Vr) (we 
again replace 167rG — > G from now on) 

^GfiVr) + ^GVt /"{Vt) + ^(2A - R) Vt /(Ft) = (89) 

with primes indicating derivatives with respect to Vr- From now on we will set the constant 
curvature density R=0; then the solutions are Bessel functions or with m = 3, 



T (90) 



MVt) = const. J3 (^aV2^Vt^ /Vi 
or 

^s{Vt) = const. Y3{(^aV2^Vt^ /V^ • (91) 

Only Jm{x) is regular as a; ^ 0, Jm{x) ~ T{m + l)~^(x/2)"^. So the only physically acceptable 
wavefunction is 

*(a, 6, . . . /) = ^{Vt) = N ^ ^^"^ ^ (92) 
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with the normahzation constant Af given by 

The latter is obtained from the wavefunction normahzation requirement 

dYrl^iVr)]^ = 1 . (94) 



f 

Jo 



Consequently the average volume of a tetrahedron is given by 

(Vt) = dVT-Vr-l *(Ft) P = "-"^ = 0.2643 . (95) 

^ ' Jo I V yi 262144\/2^/A VX ^ ' 

This last result allows us to define an average lattice spacing, by comparing it to the value for an 
equilateral tetrahedron which is Vr = (l/6\/2) ^o- then obtains for the average lattice spacing 
at strong coupling 

Zo = 1.3089 [-^ . (96) 

Note that in terms of the parameter A defined in Eq. (77) one has in all the above expressions 
VA/G = \/2A. 

The above results further show that for strong gravitational coupling, 1/G ^ 0, lattice quantum 
gravity has a finite correlation length, of the order of one lattice spacing, 

i -lo- (97) 

This last result is simply a reflection of the fact that for large G the edge lengths, and therefore 
the metric, fluctuate more or less independently in different spatial regions due to the absence of 
the curvature term. The same is true in the Euclidean lattice theory as well, in the same limit 
[23]. It is the inclusion of the curvature term that later leads to a coupling of fluctuations between 
different spatial regions. Only at the critical point in G, if one can be found, is the correlation 
length, measured in units of the fundamental lattice spacing, expected to diverge [25]. This last 
circumstance should then allow the construction of a proper lattice continuum limit, as is done in 
the Euclidean lattice theory of gravity [33] (and in many other lattice field theories as well). 



7 Perturbation Theory in the Curvature Term 

As shown in the previous section, in a number of instances it is not difficult to find the solution ^ 
of the Wheeler-DeWitt equation in the strong coupling (large G) limit, where the curvature term 
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is neglected, and only the kinetic and A terms are retained. Then the dynamics at different points 
decouples, and the wavefunction can be written as a product of relatively simple wavefunctions. It 
is then possible, at least in principle, to include the cmvaturc term as a perturbation to the zcroth 
order solution. Accordingly, the unperturturbed Wheeler-DeWitt Hamiltonian is denoted by Hq 

Ho^ - IGttG • Gij{f) + \f^) ■ 2A (98) 



dlfdli IGttG 



and the perturbation by H 



1 



1^ 'm') . (99) 

The corresponding unperturbed wavefunction is denoted by ^'o, and satisfies 

Ho^o = . (100) 
To the next order in Raleigh- Schrodinger perturbation theory one needs to solve 

(Ho + eHi) * = (101) 

where for ^ one sets as well 

^' = ^-o exple^*!} . (102) 
The sought-after first order correction *i is then given by the solution of 

iJo(*o*i) + i^i*o = 0. (103) 

Higher order corrections can then be obtained in analogous fashion. It would seem natural to search 
for a solution (here specifically in 3-1-1 dimensions) of the form 



[ <T h ) 



(104) 



with a and /? given by power series 



/T oo 

a{X,G) = ^ E ^n{GXy 



G r\ 

n=0 



^(A,G) = i^^j E/3„(GA)^ 



(105) 

The dots in Eq. (104) indicate possible higher derivative terms in the exponent of the wavefunction. 
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8 Variational Method for the Wavefunction ^ 



In this section we will describe some simple applications of the variational method for quantum 
gravity, based on the lattice Wheeler-DeWitt equation proposed earlier. The power of the varia- 
tional method is well known and appreciated in nonrelativistic quantum mechanics, atomic physics, 
and many other physically relevant applications. Its success generally rests on the ability of find- 
ing a suitable, often physically motivated, wavefunction with the lowest possible energy, thereby 
providing an approximation to both the ground state energy and the ground state wavefunction. 
In practice the wavefunction is often written as some sort of product of orbitals, dependent on a 
number of suitable parameters, which are later determined by minimization. 

Here we will write therefore an ansatz for the variational wavefunction, dependent on a number 
of free variational parameters 

= ^'[/2;a,/3,7...], (106) 

and later require that the resulting wavefunction either satisfy the Wheeler-DeWitt equation, or 
that its energy functional 

(*[^2]||-i67rG-G,,(|2)^ - [^R{l')-2\] ] Ml']) 

= ^ ^wim 

(107) 

be as close to zero as possible, |£'p = min. This procedure should then provide a useful algebraic 

relation between the variational parameters, and thus allow their determination. ^ 

Here we will consider the following correlated product variational wavefunction (in general 
dimension) 

= e-"E.^<^+^E.5,F,-h... ^ ^-1/2 -Q (e-av.^ Yl x . . . (109) 

o- h 

with variational parameters a,P, . . . real or complex. Here the J^a ^ usual volume term in d 

dimensions, and J^h ^hVh the usual Regge curvature term, in the same number of dimensions. The 
dots indicate possible additional contributions, perhaps in the form of invariant curvature squared 



®The continuum analog of the above expression would have the following general structure: 
^ ^ /d^x J[dg] <f'[g] ■ [-GA, - ^ ^{R-2X)] ■ ^[g] 

/[dff] ** [5] ■*[<;] 



(108) 



Similar energy functionals were considered some time ago by Fcynman in his variational study of Yang-Mills theory 
in 2+1 dimensions [34]. The main difference with gauge theories is that here the Hamiltonian contains two terms 
(kinetic and curvature terms) that enter with opposite signs, whereas in the gauge theory case both terms (the 
term and the (V x A)^ term) just add to each other. Feynman then argues that in the gauge theory the state of 
lowest energy corresponds necessarily to a minimum for both contributions. 
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terms. In the atomic physics literature these types of product wavefunctions are sometimes known 
as Jastrow-Slater wavefunctions [35, 36]. Note that the above wavefunction is very different from 
the ones used in minisuperspace models, as it still depends on infinitely many lattice degrees of 
freedom in the thermodynamic limit (the limit in which the number of lattice sites is taken to 
infinity) . 

The wavefunction normalization constant Z{a,j3,^ . . .) is given by 

Z = J[df]mf;a,P,...]\^ = J[dl^] expl^-2Rea^Va + 2Re/3^5hVh + ...^ (110) 

and represents the partition function for Euclidean lattice quantum gravity, but in one dimension 
less. One would expect at least Rea > to ensure convergence of the path integral; the trick we 
shall employ below is to obtain the relevant averages by analytic continuation in a and (3 of the 
corresponding averages in the Euclidean theory (for which a and P are real). Here the expression 
[c?Z^] is the usual integration measure over the edge lengths [32], a lattice version of the DeWitt 
invariant functional measure over continuum metrics [dg^,y] . The definition of Z requires that the 
functional integral in Eq. (110) actually exists, which might or might not require some suitable 
regularization: for example by the addition of curvature squared terms whose amplitude is sent to 
zero at the end of the calculation. 

Next one needs to compute the expectation value 

(*[Z2]|ii-|*[z2]) (111) 

with 

H ^ - IGttG • G,j{f) - ^ 7^ [ ^R{f) - 2A ] (112) 

2 3 

which in turn is made up of three contributions, each of which can be evaluated separately. In 
terms of explicit lattice averages, one needs the three averages, or expectation values, 

(*[a,/3,...]||-E^^(^)}l*[«'/3---]) (113) 
(*[a,/3,...]||X^y,| |*[a,/3,...]) (114) 
ma,P,...]\i^2Y,5hlh^ ma,P,...]) (115) 



with 



A,.(a) ^ Gij{f)^^. (116) 
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Note that we have summed over all lattice points by virtue of the assumed homogeneity of the 
lattice: the local average is expected to be the same as the average of the corresponding sum, 
divided by the overall number of simplices. Thus, for example, ( ^' | J^a I ^ ) = ( ^ I I ^ ) 
etc. At the same time one has, by virtue of our choice of wavefunction, 

^y,|^[a,/3,7,...]) = |*[a,/3,...]) (117) 

Y,Shlh\na,^,j,...]) = A (118) 
and also for a given simplex labeled by a 

-Ap(a)e-"^<^ = f{V,) (119) 

where / is some known function. More specifically in 2+1 dimensions one finds (here A/^ is the 
area of the relevant triangle) 

-Ai2{a)A^i = ln(n + l)Ar' (120) 

- ^"C) "^^-^ = 2531 + -TIM ('^') 



and therefore 



- Ap((7) e-°^^^ = ^ a (aAA - 2) e""^^ , (122) 



whereas in 3+1 dimensions one has (here Vt is the volume of the relevant tetrahedron) 

- (a) = 1 n (n + 6) V^'^ (123) 

7 dF VrfPF 
-A,(.)F(K.) = 1^ + ^^ (124) 



and therefore 



Ap (a) e-"^^ = 1 a (a Ft - 7) e""^^ . (125) 



In addition, in 3+1 dimensions one needs to evaluate 

- Ai2{a) e'^^h^'^^h (126) 

which is considerably more complicated. Nevertheless in 2+1 dimensions the corresponding result 
is zero, by the Gauss-Bonnet theorem. One identity can be put to use to relate one set of averages 
to another; it follows from the scaling properties of the lattice measure [dP] in d dimensions with 
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curvature coupling k = 1/SttG and unsealed eosmologieal constant Aq = A/SttG [33]. In three 
dimensions it reads 

2 Ao Vt) -k{Y, Shk) -7No = (127) 

T h 

where in the first term the sum is over all tetrahedra, and in the second term the sum is over all 
hinges (edges). The quantity A'^o is the number of sites on the lattice, the coefficient in front of 

it in general depends on the lattice coordination number, but for a cubic lattice subdivided into 
simplices it is equal to 7, since there arc seven edges within each cube (three body principals, 
three face diagonals and one body diagonal). The above sum rule can then be used by making the 
substitution Aq 2Rea and k 2Re/3. In two dimensions the analogous result reads 

2Xo{^Aa) - k{J2h)-^No = (128) 

A h 

with 2Y^f^dh = J = 47rx=const. by the Gauss-Bonnet theorem. 

Prom now on we will focus on the 2+1 case exclusively. In this case the curvature average of 
Eq. (115) is very simple 

{J^R)^47rx (129) 

where % is the Euler characteristic for the two-dimensional manifold. It will also be convenient to 
avoid a large number of factors of 167r's and make the replacement 16-jtG — > G for the rest of this 
section. Putting everything together one then finds 

_Li . + (130) 

One is not done yet, since what is needed next is an estimate for the average area of a triangle, 
^A- This quantity is given, for a general measure over edges in two dimensions of the form 

^d/2.^T(^A)^by 

(^a) = , (131) 



again with the requirement Rea > for the average to exist. It will be convenient to just set in 

2 
3 



the following Aa = {Aa) = ctq/q; with ao = (1 + I cr)/2. One then obtains, finally, the relatively 



simple result 



E{a) _ (70 - 2 2Aao 1_ _ Attx 

GNt ~ 4 ■ " + G2 ■ a G^T " ^ ' 



It would seem that, in order to avoid a potential instability, it might be safer to choose ao > 2. 
The roots of this equation (corresponding to the requirement ( ^' | i7 | ^' ) = ) are given by 
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with 

A = QA-k'^x^ - 8G^N^Xao{ao - 2) , (134) 

so that A is zero for 

G = Ge = ± '^"^ , . (135) 

Here we select, on physical grounds, the positive root. When A = (or G = Gc) the two complex 
roots become real, or vice-versa, with 

ao(Gc) = EE^^ > if X > . (136) 

Thus for strong couphng (large G > Gc) a is almost purely imaginary 

whereas for weak coupling (small G < Gc) the two roots become 

ao ^ ^^ + OiG^) 

"° ^ G2iVT(ao - 2) - + 0{G). (138) 

Note that an identical set of results would have been obtained if one had computed |£'(a)p for 
complex alpha, and looked for minima. This is the quantity displayed in Figures 6 and 7. 

Next we come to a brief discussion of the results. One interpretation is that the variational 
method, using the proposed correlated product wavefunction in 2+1 dimensions, suggests the pres- 
ence of a phase transition for pure gravity in G, located at the critical point G = Gc- This picture 
found here would then be in accordance with the result found in the Euclidean lattice theory in 
Ref. [31], which also gave a phase transition in three-dimensional gravity between a smooth phase 
(for G > Gc) and a branched polymer phase (for G < Gc)- A similar transition was found on the 
lattice in four dimensions as well [23]. Finally, the presence of a phase transition is also inferred 
from continuum calculations for pure gravity in e = d — 2 > 2, although the latter does not give a 
clear indication on which phase is physical; nevertheless simple renormalization group arguments 
suggest that the weak coupling phase describes gravitational screening, while the strong coupling 
phase implies gravitational anti-screening. This last expansion then gives a critical point for pure 
gravity in 2+1 dimensions at Gc = ^(^^ — 2) + ^Q^d — 2)^ + . . ., or Gc ~ 0.024 in units of the 
cutoff [44, 45, 46]. The Euclidean lattice calculation quoted earlier gives, in the same dimensions, 
Gc ~ 0.355. Note that the numerical magnitude of the critical point G in lattice units, contrary to 
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the critical exponents, is not expected to be universal, and thus cannot be compared directly be- 
tween formulations utilizing different ultraviolet regulators. We shall not enter here into some of the 
known peculiarities of three-dimensional gravity, including the absence of perturbative transverse- 
traceless radiation modes, and the absence of a sensible Newtonian limit; a recent discussion of 
these and related issues can be found for example in [25], and further references cited therein. 

In 3-1-1 dimensions the variational calculation is quite a bit more complex, since the integrated 
curvature term is no longer a constant. In the small curvature limit and for small variational 
parameter /3 we have obtained the following expansion for the variational parameter a 



Here cq is a real constant whose value we have not been able to determine yet. The two roots are 
found to become degenerate and real for 



which is again suggestive of a phase transition at Gc in 3-1-1 dimensions, as found previously in 
the Euclidean theory in four dimensions [23, 33]. More detailed calculations in the 3-1-1 case are in 
progress, and will be presented elsewhere [47]. 

We conclude this section by observing that our results suggest a rather intriguing relationship 
between the ground state wave functional of quantum gravity in n + 1 dimensions, and averages 
computed within the Euclidean Feynman path integral formulation in n dimensions, i.e. in one 
dimension less. Moreover, since the variational calculations presented here rely on what could 
be regarded as an improved mean field calculation, they are expected to become more accurate 
in higher dimensions, where the number of neighbors to each lattice point (or simplex) increases 
rapidly. 




(139) 




(140) 
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Figure 6. Energy surface |i?(a)p in 2+1 dimensions at strong coupling, G ^ Gc in the {Rea, Ima) plane. 
Note the presence of two almost purely imaginary, complex conjugate roots. The specific values used here 
are x = 2, Nt = 10, ao = 3 and A = 1. 




Figure 7. Energy surface |i?(Q:)p in 2+1 dimensions for weak coupling, G ^ Gc- In this case both roots are 
along the real a axis. 
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9 Weak field expansion 



In this section we will discuss briefly the weak field expansion for the proposed lattice Wheeler- 
DeWitt equation. The purpose here is to show how the weak field expansion is performed, and how 
results analogous to the continuum ones are obtain for sufficiently smooth manifolds. Such results 
would be of relevance to the weak coupling (small G) expansion, and to an application of the WKB 
method on the lattice, for example. More generally a clear connection to the continuum theory, and 
thTis between lattice and continuum operators, is desirable, if not essential, in order to understand 
the meaning of physical gravitational averages, such as average curvature etc. First we note here 
that the lattice kinetic term (the one involving Gij) has the correct continuum limit, essentially 
by construction. On the other hand the curvature term appearing in the discrete Wheeler-DeWitt 
equation in 3+1 dimension is nothing but the integrand in the Regge expression for the Einstein- 
Hilbert action in three dimensions, 

Ie = -k Yl ^hSh ■ (141) 

edges h 

The expansion of this action around fiat space was already considered in some detail in Ref. [31], 
and shown to agree with the weak field expansion in the continuum. Here we provide a very short 
summary of the methods and results of this work. Following Ref. [20], one takes as background 
space a network of unit cubes divided into tetrahedra by drawing in parallel sets of face and body 
diagonals, as shown in Figure 8. With this choice, there are 2*^ — 1 = 7 edges per lattice point 
emanating in the positive lattice directions: three body principals, three face diagonals and one 
body diagonal, giving a total of seven components per lattice point. 
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Figure 8. A cube divided into simplices. 

It is convenient to use a binary notation for edges, so that the edge index corresponds to the 
lattice direction of the edge, expressed as a binary number 

(0,0,1) ^ 1 (0,1,1)^3 (1,1,1)^7 
(0,1,0) ^ 2 (1,0,1)^5 

(1,0,0) ^ 4 (1,1,0)^6 (142) 

The edge lengths are then allowed to fluctuate around their flat space values Zj = Z°(l + e^), and 
the second variation of the action is expressed as a quadratic form in e 

^2 J ^ J2 e("*) ^ M("*'") e^'*), (143) 

mn 

where n,m label the sites on the lattice, and is some Hermitian matrix. The general aim 

is then to show that the above quadratic form is equivalent to the expansion of the continuum 
Einstein-Hilbert action to quadratic order in the metric fluctuations. The infinite-dimensional 
matrix M^"*'") is best studied by going to momentum space; one assumes that the fluctuation ej 
at the point j steps from the origin in one coordinate direction, k steps in another coordinate 
direction, and I steps in the third coordinate direction, is related to the corresponding fluctuation 
€i at the origin by 

(j+k+l) _ , J- 1 M (-,AA -\ 

with coi = e^^^. In the smooth limit, ooi = 1 + iki + 0{kf), the lattice action and the continuum 
action are then expected to agree. Note also that it is convenient here to set the lattice spacing in 
the three principal directions equal to one; it can always be restored at the end by using dimensional 
arguments. 

It is desirable to express the lattice action in terms of variables which are closer to the continuum 
ones, such as /i^jy or h^i, = h^y — '^rj^uhw. Up to terms that involve derivatives of the metric 
(and which reflect the ambiguity of where precisely on the lattice the continuum metric should be 
defined), this relationship can be obtained by considering one tetrahedron, and using the expression 
for the invariant line element ds^ = g^ydx^dx'^ with g^y = rj^y + h^y, where rj^y is the diagonal 
flat metric. Inserting Zj = /° (1 + e^), with Z° = 1, \/2, \/3 for the body principal {i = 1, 2,4), face 
diagonal {i = 3,5,6), and body diagonal {i = 7), respectively, one obtains 

(l + ei)2 = 1 + hu 
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(1 + 62)' 


= 1 + /i22 




(1 + 64)' 


= 1 + /l33 




(1+^3)' 


= 1 + ^(^11 


+ h22) + hi2 


(l + es)' 


= ^ + \{hii 


+ /i33) + ^13 


(l + ee)' 


= l + \{h22 


+ hz2,)+ h23 


(1 + e^f 


= l + \{hii 


2 

+ h22 + /i33) + ^(/ll2 + /i23 + /il3 



(145) 



(note that we use the binary notation for edges, but maintain the usual index notation for the field 
h^iy). The above relationship can then be inverted to give the e's in term of the /I's . Note that 
there are seven variables, but only six h^j^^s (in general in d dimensions we have 2*^ — 1 variables 
and d{d + l)/2 /i^jy's, which leads to a number of redundant lattice variables for d > 2). 

Thus, to lowest order in h^^, one can perform a field rotation on the lattice in order to go from 
the Ci variables to the /ly^^i/'s (or huvs), 



(146) 



with 



e = Uih h = U2h , 



(147) 



and so 



e = Vh V = U1U2 . 
Here V and Ui are 7x6 matrices, while C/2 is a 6 x 6 matrix, 
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1/ 



(148) 



(149) 



(150) 



The above rotation is an essential step in transforming the lattice action into a form that looks 
like the continuum action, to quadratic order in the weak fields. For the Regge-Einstein term the 
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matrix M^^ describing the small fluctuations around flat space is given by 



(M,)i,i = 

(M^)i,2 = 

(M^)l,4 = 
(M^)l,6 = 
(M^)l,7 = 
(M^)4,4 = 
(M^)4,5 = 
(M^)4,7 = 



-2 

—UJliU4 — U}204 

2 + 2^2 
2coi + 2ci;2Ct)4 
— 80)2 — 80)4 
-8 

— 4a;2 — 4a)4 
6 + 6u4 
-18 



(151) 



with the remaining matrix elements obtained by permuting the appropriate indices. Due to its 
structure, which is of the form 

M^=(^f ^.V (152) 



-18 



where is a 6 x 6 matrix, a rotation can be done which completely decouples the fluctuations in 
£7, 



with 



^6 + 4^bb^ 



i6 I 






-18 



h 
1 



(153) 

(154) 



One then finds the first important result, namely that the small fluctuation matrix has three 
zero eigenvalues, corresponding to the translational zero modes for M^^ 
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\ i(l - a;ia;2W4) 
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i(l - a;2a;4) 






1 — bJ\ 





^(1 



2(1 - a;2a;4) 

i(l - U;iU2UJ4:)/ 




(155) 



^(1 — UJ1UJ2UJ4) 

where xi,X2,X2 are three arbitrary parameters. The remaining eigenvalues are —18 (once) and 
0{k^) (three times). Notice that one mode, corresponding to the fluctuations in the body diagonal 
er, completely decouples. The next step is to transform the lattice weak field action into a form 
similar (in fact identical) to the continuum form. One further rotation by the (6 x 6) matrix T^, 
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defined by 



\ 

gives the new small fluctuation matrix 
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M' T 



(157) 



This last transformation is equivalent to a change to trace-reversed metric variables. Finally one 
defines the gauge fixed matrix 

= L^- \ciC^, (158) 

where C^^ is introduced in order to give the lattice analog of the harmonic gauge fixing term, with 
/5(-l+a;i) 1-^1 1-u;, 6 (l - 4) 6 (l - ^) \ 

C^ = l 5(-l + u;2) 1-UJ2 6(1-^) 6(1-^) .(159) 

V l-c^4 5(-l + u;4) 6(1-^) 6(l-^)y 

Then the form of L^^ is precisely equivalent to the corresponding continuum expression, in trace- 
reversed variables and in the harmonic gauge [31]. The seemingly complex combined and 
rotations just correspond to a rotation, from the original lattice edge fluctuation variables (e) to 
the trace reversed metric variables (h). 

Perhaps the most important aspect of the above proof of convergence of the lattice curvature 
term, and more generally of the whole lattice Wheeler-DeWitt equation, towards the corresponding 
continuum expression is in its relevance to the weak field limit, to a perturbative expansion in G, 
and to a WKB expansion of the wavefunction. The latter are all issues that have already earned 
some consideration in the continuum formulation [4, 6, 7, 10]. The present work suggests that most 
of those continuum results will remain valid, as long as they are derived in the context of the stated 
approximations. In particular, there is no reason to expect the lattice semiclassical wavefunction 
to have a different form (apart from the use of different variables, whose correspondence has been 
detailed in this section) compared to the continuum one [7] . 



10 Conclusions 



In this paper we have presented a lattice version of the Wheeler-DeWitt equation of quantum 
gravity. The present 3-1-1 approach is based on the canonical, and therefore Lorentzian, formulation 
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of quantum gravity, and should therefore be regarded as complementary to the four-dimensional 
Euclidean lattice version of the same theory discussed earlier in other papers. The equations are 
explicit enough to allow a number of potentially useful practical calculations, such as the strong 
coupling expansion, mean field theory, and the variational method. In the preceding sections we 
have outlined a number of specific calculations to illustrate the mechanics of the lattice theory, 
and the likely physical interpretation of the results. Because of its reliance on different set of 
nonperturbative approximation methods, the formulation presented here should be useful in viewing 
the older Euclidean lattice results from a very different perspective; in a number of instances we 
have been able to show the physical similarities between the two types of results. 

Nevertheless the phenomenal complexity of the original continuum theory, and of the Euclidean 
lattice approach, with all its issues of, for example, pcrturbative nonrenormalizability in four di- 
mensions, has not gone away; it just got reformulated in a rather different language involving a 
Schrodinger -like equation, wavefunctionals, operators and states. The hope is that the present 
approach will allow the use of a different set of approximation methods, and numerical algorithms, 
to explore what in some instances are largely known issues, but now from an entirely different 
perspective. Among the problems one might want to consider we list: the description of invariant 
correlation functions [37, 38] , the behavior of the fundamental gravitational correlation length ^ as 
a function of the coupling G, the approach to the lattice continuum limit at Gc, estimates for the 
critical exponents in the vicinity of the fixed point, and the large-scale behavior of the gravitational 
Wilson loop [39]. Note that what is sometimes referred to as the problem of time does not neces- 
sarily affect the above issues, which in our opinion can be settled by looking exclusively at certain 
types of invariant correlations along the spatial directions only. These Green's functions should 
then provide adequate information about the nature of correlations in the full theory, without ever 
having to make reference to a time variable. 
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Appendix 



A The matrix M in 3+1 dimensions 

The matrix of coefficients of the second partial derivative operators in 3+1 dimensions is given 
by ^^7, where M is a symmetric 6x6 matrix with entries as follows: 

Mil = 

Mi2 = a^ + b^ + c^ - 2ab - 2b 

Mi3 = a^ + cf + 6"^ - 2ae - 2de 

Mi4 = a^ + b^ + - 2ac - 2bc 

Mi5 = + (f + - 2ad - 2de 

Mi6 = 6^ + + + - 2af + 2bd + 2ce - 2bc - 2be - 2cd - 2de 

M22 = 2c2 

M23 = (? + d^ + f-2cf-2df 

M24 = + b^ + - 2ab - 2ac 

M25 = + 6^ + + /2 + 2a f + 2bd - 2ce - 2ab - 2ad - 2b f - 2df 

M26 = + d^ + f -2cd-2df 

M33 = 2d^ 

M34 = + + e+/^ + 2af - 2bd + 2ce - 2ac - 2ae - 2c/ - 2e/ 

M35 = + d^ + - 2ad - 2ae 

M36 = + d^ + f -2cd-2cf 

M44 = 2b^ 

M45 = b'' + e^ + f-2bf-2ef 

M46 = b^ + + f - 2be-2ef 

M55 = 

M56 = b^ + + f --2be-2bf 

Mm = 2/2. (1) 
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B Lattice Hamiltonian for Gauge Theories 



It is of interest to see how the Hamiltonian approach has fared for ordinary SU{N) gauge 
theories, whose nontrivial infrared properties (confinement, chiral symmetry breaking) cannot be 
seen to any order in perturbation theory, and require therefore some sort of nonperturbative ap- 
proach, based for example on the strong coupling expansion. In the continuum one starts from the 
Yang-Mills action 



I = -^^jd'xF^^F^^'^ (1) 

with field strength 

F;, = d,Al - d,Al + A^Al (2) 

and gauge fields A^ with (a = 1 . . . A'"^ — 1), where the quantities fahc are the structure constants 
of the Lie group, such that the generators satisfy [TajT;,] = ifabcTc- 

A lattice regularized form of the gauge action in Eq. (1) was given in Ref. [40], see also [41]. 
The theory is defined on a d-dimensional hyper-cubic lattice with lattice spacing a, vertices labeled 
by an index n and directions by /x. The group elements = eyip{iag A^^Ta) are defined in 
the fundamental representation, and reside on the links of the lattice. The pure gauge Euclidean 
action involves a sum of traces of path-ordered products [with C/_^(n + v) = UJ^{n)] of unitary 
Un{n) matrices around an elementary square loop (" plaquettes" , here denoted by □), 

I[U] = ^tr [t/C/C/tf/t + h.c] . (3) 

9 □ 

The action is locally gauge invariant with respect to the change 

C/^ (n) ^ V\n)U^{n)V{n + i^) , (4) 

where V is an arbitrary SU{N) matrix defined on the lattice sites. 
The next step is to define the path integral as 

Z{g'') = jidUn] exp (-/[[/]) , (5) 

where [dUn] is the Haar measure over the group SU{N), one copy for each lattice link variable 
U. A lattice regularized Hamiltonian can then be defined on a purely spatial lattice, by taking the 
zero lattice spacing limit in the time direction [42, 43]. Local gauge invariance further allows one 
to set all the link variables in the time direction to unity, C/„o = 1, or A^q = in this lattice version 
of the temporal gauge. The U variables can now be eliminated by introducing generators of local 
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rotations Ef{n), defined on the links (with spatial directions labeled by i,j = 1, 2, 3) and satisfying 
the commutation relations 

[Ef (n), Uj{ra)] = Ui{n) dij (5„™ , (6) 
along with the SU {N) generator algebra commutation relation 

'E!l{n), 4(m)] = i f^'^ Efin) 5^^ 5^^ . (7) 

This finally gives for the Hamiltonian of Wilson's lattice gauge theory [42] . 

H = E'^E" - y — ^ tr \UUU^U^ + h.c] (8) 

links Q ^ 

The first term in Eq. (8) is the lattice analog of the electric field term E^, while the second term is 
a lattice discretized version, involving lattice finite differences, of the magnetic field (V x A)^ term. 
In this picture the analog of Gauss's law is a constraint which needs to be enforced on physical 
states at each spatial site n 

yEt{n)\^)=0 (9) 

i=l 

In general, and irrespective of the symmetry group chosen, the ground state in the strong coupling 
g'^ oo limit has all the SU{N) rotators in their ground state. In this limit the Hamiltonian has 
the simple form 

^0 = E EfEf . (10) 

links 

Then the vacuum is a state for which each link is in a color singlet state 

e: io) = . (11) 

The lowest order excitation of the vacuum is a state with one unit of chromo-electric field on each 
link of an elementary lattice square, and energy 

g^ N'^-l 

En = . (12) 

Raleigh-Schrodinger perturbation theory can then be used to compute corrections to arbitrarily high 
order in l/g^. But ultimately one is interested in the limit g^ — > 0, corresponding to the ultraviolet 
asymptotic freedom fixed point of the non-Abclian gauge theory, and thus to the lattice continuum 
limit a —7- 0. Thus in order to recover the original theory's continuum limit, one needs to examine 
a limit where the mass gap in units of the lattice spacing goes to zero, am{g) — >■ 0. This limit 
then corresponds to an infinite correlation length in lattice units; the zero lattice spacing limit so 
described is a crucial step in fully recovering desirable properties (rotational or Lorentz invariance, 
asymptotic freedom, massless perturbative gluon excitations etc.) of the original continuum theory. 
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